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Abstract. 

We establish a matrix generalization of the ultradiscrete fourth Painleve equation (ud-Prv). 
Well-defined multicomponent systems that permit ultradiscretization are obtained using an 
approach that relies on a group defined by constraints imposed by the requirement of a consistent 
evolution of the systems. The ultradiscrete limit of these systems yields coupled multicomponent 
ultradiscrete systems that generalize ud-Piv- The dynamics, irreducibility, and integrability of 
the matrix valued ultradiscrete systems are studied. 



AMS classification scheme numbers: 39A13, 33C70, 37J35, 16Y60 
1. Introduction 

Discrete Painleve equations are difference equation analogs of classical Painleve equations [32] 
and have been extensively studied recently (see the review article [6]). The ultradiscrete Painleve 
equations are discrete equations considered to be extended cellular automata (they may also 
be considered as piecewise linear systems) that are derived by applying the ultradiscretization 
process [21] to discrete Painleve equations. This process has been accepted as one that preserves 
integrability [7J. Particular indicators of integrability in the ultradiscrete setting include the 
existence of a Lax pair [19], an analog of singularity confinement [9], and special solutions [16J. 

All of the preceding examples arising from ultradiscretization are one-component (that 
is, scalar) systems. Generalizations of integrable systems to associative algebras have been 
considered for many years (see [TSj, and references therein). However, the general methods and 
results previously obtained are inapplicable in the ultradiscrete setting, due to the requirement of 
a subtraction free setting. We present for the first time a matrix generalization of an ultradiscrete 
system. 

The constraints related to the subtraction free setting and consistent evolution are studied 
in a group theoretic approach, in which one may also describe the nature of the irreducible 
subsystems. As an application of this method, we introduce a matrix version of the ud-Prv of 
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[TT] which is derived by applying the ultradiscretization procedure to q-Prv, 

f ( n +\ — „ „ f (+\ 1 + ffl 2/2(f)+a2ao/2(f)/o(t) 

JoW) - «oai/i^J 1+ao/oW+aoai/oW/i(t) 



f ( n -l-\ — n n f (+\ l+ao/oW+ao«l/oW/l(<) /i i\ 

h(qt) - aia 2 / 2 (t) 1+ai/i(i)+aia2/i(t)/2(t) (1.1) 

f ( n +\ — n n f (+\ l+ai/ift)+ma2/i(f)/2(t) 

With the explicit form of the matrices derived, the new systems can be considered as coupled 
multicomponent generalizations. It should be stressed that the approach of this paper gives all 
possible ultradiscretizable matrix valued versions of (11.11) . 

The reason for choosing g-Prv is that it has already been thoroughly and expertly 
investigated in the scalar case (i.e., when /j and cij are scalar) [II] . In [TTj . g-Ppv was 
shown to admit the action of the affine Weyl group of type A 2 as a group of Backlund 
transformations, to have classical solutions expressible in terms of g-Hermite- Weber functions, 
to have rational solutions, and its connection with the classification of Sakai [23] was also 
investigated. Furthermore, the ultradiscrete limit was taken in [TT], and was shown to also 
admit affine Weyl group representations. As g-Prv is such a rich system, and has already been 
well-studied, this makes it a perfect system for the application of our approach of ultradiscrete 
matrix generalization. 

Before turning to the derivation of matrix ud-Piv, ultradiscretization should be introduced 
in more detail, so that the reason for certain constraints given later will be clear. 

The process is a way of bringing a rational expression, /, in variables (or parameters) 
ax, . . . , a n to a new expression, F, in new ultradiscrete variables A\, . . . , A n , that are related to 
the old variables via the relation aj = e Ai ^ e and limiting process 

F(A u ...,A n ) =\ime\ogf(a u ...,a n ). (1.2) 

e— >0 

In general it is sufficient to make the following correspondences between binary operations 

a + b — > max(yl, B) 
ab -> A + B (1.3) 
a/b A-B. 

This process is a way in which we may take an integrable mapping over the positive real 
numbers R + to an integrable mapping over the max-plus semiring [5]. The requirement that the 
pre-ultradiscrete equations are subtraction free expressions of a definite sign is a more stringent 
restraint in the matrix setting than the one-component setting, and it is this requirement which 
motivates the particular form of our matrix system. 

The outline of this paper is as follows. In section [21 a g-Prv is derived in the noncommutative 
setting, where the dependent variables take their values in an associative algebra. In section [3] 
conditions on the matrix forms of the dependent variables and parameters of g-Prv are derived 
such that it has a well-defined evolution and is ultradiscretizable. The group theoretic approach 
is adopted to describe the constraints on the system. In section 0] the ultradiscrete version of this 
system is derived, and some of the rich phenomenology of the derived matrix valued ultradiscrete 
Piv is displayed and analyzed in section [5] 
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2. Symmetric g-Prv on an associative algebra 



In this section it is shown that the symmetric g-Prv of [H] can be derived from a Lax formalism 
in the noncommutative setting, where the dependent variables take values in an a priori 
arbitrary associative algebra, A, with unit I over a field K (when we turn to ultradiscretization, 
the requirement of a field will be modified, but not in such a way as to affect the derivation 
from a Lax pair). This puts the present work in the context of other recent work on integrable 
systems such as [H] and [15] where the structure of integrable ODEs and PDEs (respectively) 
was extended to the domain of associative algebras, and [1] where Painleve equations were 
defined on an associative algebra (see also [15]). This trend has also been present in work on 
discrete integrable systems, such as [2] where the higher dimensional consistency (consistency 
around a cube) property was investigated for integrable partial difference equations defined on 
an associative algebra, and jl] where an initial value problem on the lattice KdV with dependent 
variables taking values in an associative algebra was studied, leading to exact solutions. 

The auxiliary (spectral) parameter x, time variable t and constant q belong to the field IK. 
The dependent variables {fi} G A, system parameters G A {i G {0, 1,2}), and we define 



r l :=/ + ^ /i + & 3 6 3 +i/i/i+i; (21) 

up to an arbitrary ordering of the and fj factors. (It will be shown that the ordering of these 
factors within Tj is of no consequence for either the integrability of the system or the existence 
of a well defined evolution in the ultradiscrete limit.) The invertibility of these expressions is 
assumed, that is T^ 1 G A. 

We derive the system from a linear problem to settle other ordering issues in the 
noncommutative setting. The g-type Lax formalism is given by 



4>(qx,t) = L(x,t)(f)(x,t) , (f>(x,qt) = M(x,t)(j)(x,t) 



(2.2) 



where 



L(x,t) 



M(x,t) 



hfat a 
[l + x 2 )I 






bifit'l 
{l + x 2 )I 

\ 



\ 



J 



(2.3a) 



(2.36) 



;i + x 2 )/ 
o 

/ b 2 2 T 2 

& - 2 r 
\&r 2 ri j 

The ultradiscrete version of this linear problem (for the usual commutative case) originally 
appeared in [TO] . 

The compatibility condition for this linear problem reads 

M(qx, t)L(x, t) = L(x, qt)M(x, t), (2.4) 



and leads to 



61/1 
hf 2 



q ^ bf r 2 h h r 



-1 1,2 



<? 2/3 & - 2 ro^rr 1 ^ 

q^bfT^hT-^h: 



(2.5) 
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where the overline denotes a time- update and 6j = b^. 

Following [11], we show a product of the dependent variables can be regarded as the 
independent variable. With {f^ 1 } G A, {b^ 1 } G A (i.e., we are working with a skew field) and 
specifying that the product &0/0&1/1&2/2 is proportional to /, it is seen that = QC 2 I 

where c G K and c = qc. Without loss of generality we set c = t. From now on 

60/0&1/1&2/2 = qt 2 I (2.6) 

will be imposed (so the algebra generated by all three {f{\ and I is not free). The invertibility 
of the algebra elements {f{\ and {b{\ is a consequence of the explicit matrix representation of 
these objects for the well-defined matrix systems studied in the next sections. 
With the restriction 

blb\b\ = ql (2.7) 

imposed, the map (12.51) is a noncommuting generalization of g-Prv- If we specify that bi and fi be 
matrix valued, the only requirement for a consistent evolution is that the Tj are invertible. This 
however is too general a system to be ultradiscretized, since in general we require the inverse to 
be subtraction free. 

If all variables commute, then after the change of variables a« := 6f the map reduces to 
g-Piv, (11 .ip . as presented in [TT] . 



3. Ultradiscretizable matrix structure 

The conditions (12.61) and (12. 7p can be used in conjunction with (12.51) to define constraints that 
lead to a consistent evolution on A as a free algebra with two constant (say b\ and 62) and two 
variable (say f\ and $2) generators. Regarding these as n x n (or even infinite dimensional) 
matrices leads to multicomponent systems. However, the aim of the present work is to derive 
matrix (or multicomponent) ultradiscrete systems, and hence, as we require the expressions to 
be subtraction free, we have considerably less freedom than this general setting. 

Due to this restriction, we restrict A to be the group of invertible non-negative matrices, 
that is we set 

A = S n x K n (3.1) 

where S n is the symmetric group and K will further be restricted to be M + in models where we 
wish to perform ultradiscretization. For our purposes S n is realized as n x n matrices of the 
form ^io-(j) for a G S n . (This group decomposition result can been seen in [3].) We define the 
homomorphism 71 : A — > A/KJ 1 = S n to be the homomorphism obtained as a result of the above 
semidirect product. This allows us to more easily deduce the form of the matrices {fi}, {h}, 
that give a well-defined evolution. 

Since A is a semidirect product, the elements bi and fi can be uniquely written in the form 



bi Pi Si 

fiit) = Fi(t)zi, 



(3.2) 
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where 7r(6j) = G S^, 7r(/j(i)) = ^ G S" n , and $ and F j(t) are diagonal matrices containing the 
n components of 6j and fi(t) respectively (we leave the matrix representation implicit). 

We now derive further restrictions on {si} and {z{\ such that the evolution is consistent, 
and all terms in the map (such as the Tj) remain in A, (13. ip . 

Consider the following form of T,, 

r, : / • !>!.[, ■ W.^L-lI). (3.3) 

As 7r(J) = /, Tc(Ti) = I and this implies 

s\zi = I % G {0,1, 2}. (3.4) 

This is the only condition that arises from the requirement that Tj G A, where A is given by 
(13. ip . It is immediately seen that condition (13.41) is independent of the ordering of the b\fi 
term in IV There are 24 possible orderings of b^bf +1 f i+1 fi (we do not consider the possibility of 
splitting up the 6, factors, as bf =: a, is the parameter in the commutative case [H]). Of these 24 
possibilities, 8 also require the commutativity of and s| (equivalently Zj and z i+ i). It is shown 
in Appendix A that these additional commutativity relations do not change the restrictions on 
{si} and {zi}. (That is, commutativity of sf and is a consequence of the full set of relations.) 

Requiring the preservation of (13. 2p as the variables evolve, the projection of (12. 5p onto S n , 
with (13. 4p . gives 

sf = 3^1, i€{0,l,2}. (3.5) 
The projection of the constraints (I2.6P and (12 .7p onto S n , with (13 .4p . gives 

sjsjsl = I (3.6) 

and 

444 = i (3.7) 

respectively. 

Therefore, to give a consistent evolution that permits ultradiscretization, {s{\ are 
homomorphic images of the group generators of 

ri I l 4 224 224 22222 i 333 i\ /oo\ 

G= {9o,9u92\9o =9x92, 9i = 9 2 9 i 9 2 = 9o9i, 9 2 9i9 = !> 9o9i9 2 = l ) ( 3 -8) 
in S n ; {z^ are given by (13.41) . The group G has order 108. The order of the generators of G is 



shown to be 18 in Appendix A 



4. Ultradiscretization 



We now consider the ultradiscretization of the matrix valued systems derived in the previous 
section. The components of the ultradiscretized systems belong to the max-plus semiring, S, 
which is the set RU{— oo} adjoined with the binary operations of max and + (often called tropical 
addition and tropical multiplication). To map the pre- ultradiscrete expression to the max-plus 
semiring, we may simply make the correspondences (11.31) on the level of the components. (So — oo 
becomes the additive identity and becomes the multiplicative identity.) By ultradiscretizing 
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max(ajj, bif) 
max(a tfc + b kj ) 

k 



matrix operations, we arrive at the following definitions of matrix operations over S. If A = (a^) 
and B = (fey), then following [19], we define tropical matrix addition and multiplication, © and 
©, by the equations 

(A © B) %3 
(A © B)ij 

along with a scalar operation given by 
(A © A)ij := (A + otj) 

for all A G S. In the ultradiscrete limit is mapped to — oo, and 1 is mapped to 0; hence the 
identity matrix, /, is the matrix with 0s along the diagonal and — oo in every other entry. In the 
same way it is clear what happens to matrix realizations of members of S n in the ultradiscrete 
limit. 

An ultradiscretized member of the group A, ( 13. ip . has a decomposition of the form 
D = A ©T 

(cf. equation (13.21) ) where A has — oo for all off-diagonal entries and T is an ultradiscretization 
of an element of S n . Its inverse is given by 

zr 1 = T~ l © A -1 , 

where (A _1 )jj = —(A) u and all off-diagonal entries are — oo. 

As well as the matrix map, the correspondence also allows us to easily write the Lax pair 
over the semialgebra. 

( max(0, 2X) © I B © F © -oo \ 

— oo 



C(X,T) 



M(X,T) 



\ 



Bo 



F 2 ©-^ 



max(0,2X)©7 B x © F x © -f- 



-oo 



-oo 
-oo 



B. 



2 « 1 2 

— OO 

— OO 



-OO 



max(0,2X) 
\ 



/ 



Bo 2 



To 



l B x © Ti — oo — oo 

Where the ultradiscretization of Tj, as given in ( 12.11) . is the matrix 



? 3 



F i+ i) • 



r, = / © (Bf © Ff) © (Bf © Bl 

The compatibility condition reads 

M{X + Q, T) © C{X, T) = C(X, T + Q)®M(X,T), 
and gives the ultradiscrete equation over an associative S'-algebra 



B ®F 
B\ © F\ 
B 2 ®F 2 



\Q®B 2 



r 2 © B x © Fx © Tr 



3Q « ' A :i 



5 n 



)-2 



r ©5 2 

Ti © 5 



F 2 



rr 1 
rr 1 



#0) 

2 

2' 



The ultradiscrete version of the restrictions (12.61) and (12.71) are 

B © F © B x © Fx © 5 2 © F 2 = (Q + 2T) © J 



(4.1a) 

(4.16) 

(4.2) 
(4.3) 

(4.4) 

(4.5) 
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Rank 


On 

yu 


fih 

yi 


ill 


1 


1 


1 


1 


2 


fl 2) 


fl 2) 


1 


3 


fl 2 3) 


fl 2 3) 


fl 2 3) 




(1,2,3) 


(1,3,2) 


1 


4 


(1,2)(3,4) 


(1,3)(2,4) 


(M)(2,3) 


6 


(1,2,3,4,5,6) 


(1,6,5,4,3,2) 


1 




(1,2)(3,4)(5,6) 


(1,3,6)(2,4,5) 


(1,5,3,2,6,4) 




(1,2,3)(4,5,6) 


(1,4,3,6,2,5) 


(1,4,3,6,2,5) 



Table 1. Lowest rank cases of homomorphic images of the generators of G in S, 



and 

B 3 ® E 3 ® E 3 = Q <g> /. (4.6) 

(Of course, it would have been equally legitimate to apply the correspondence on the level of the 
map (12.51) without starting from a derivation from the ultradiscretized Lax pair.) 

It is easily seen that if 2Q/3, the parameter T, and all components of the map belong to Z 
then at all time-steps all components (not formally equal to — oo) belong to Z. It is this property 
which motivates the term 'extended cellular automata'. 

5. Phenomenology 

As mentioned in the above discussion, we are required to find homomorphic images of the group 
G in S n . To do this, we use the computer algebra package Magma. The homomorphic images 
of G in S n give rise to reducible and irreducible subgroups, which in turn translate to reducible 
and irreducible matrix valued systems. By definition, the reducible systems are decomposable 
into irreducible systems, and hence we restrict our attention to the irreducible cases. 

We may use any homomorphism to induce a group action of G onto a set of n objects. In 
this manner, we may state by the orbit stabilizer theorem that the size of any orbit of G must 
divide the order of the group. Since the group has order 108, this implies the irreducible images 
of G be of sizes that divide 108. In terms of matrix valued systems, the implication is that any 
irreducible matrix valued systems are of sizes that divide 108. 

The lowest rank cases of the homomorphic images of the generators of G in S n are given 
in table [1] using the standard cycle notation for the symmetric group. The rank 1 case is well 
understood [IT] ; hence we turn to the rank 2 case. For the examples presented here, we restrict 
our attention to the ordering within the {T,}, (14.21) . 

Typical behavior of the rank 2 map is shown in figure [TJ The initial conditions and parameter 
values in this case are 

*-(7-L) Mt-°~) 
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where B2 and F2 are determined by the constraints, and Q = 1. For most initial conditions and 
parameter values, the behavior has a similar level of visual complexity. 




(a) Fo: both components (b) Fi; both components (c) F2: both components 

Figure 1. Generic behavior of the rank 2 case. Component values are plotted against time step 
values. 



It is a hallmark of the integrability of Painleve systems that they possess special solutions 
such as rational and hypergeometric functions [2T] . A remarkable discovery of our numerical 
investigations is that (14.41) displays special solution type behavior. These solutions only occur 
for specific parameter values and initial conditions. One example of this comes at a surprisingly 
close set of parameters and initial conditions to those displayed by figure HJ By setting the 
parameters to be 



-00 





Bi 



—00 

3 
5 




with the same set of initial conditions, the behavior coalesces down to the much simpler form 
shown in figure [21 




(a) Fo: both components (b) F\: both components 

Figure 2. Some special behavior of the rank 2 case. 



(c) F2: both components 



The graphs of the single components in figure [2] strongly resemble the recently discovered 
ultradiscrete hypergeometric functions of [16J. This implies that the special solution behavior 
shown here may be parameterized by a higher-dimensional generalization of the ultradiscrete 
hypergeometric functions of [16]. We discuss this possibility further in section [6j Behavior 
resembling rational solutions has also been observed in our computational investigations. 
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The typical behavior of the rank 3 map is shown in figure [3J The initial conditions and 
parameter values are 
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— oo 
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— oo 
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— OO 


-5 


— oo 
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— oo 


— oo 
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— OO 


— oo 


1 / 



where the coupling comes from the forms of the parameters. 

I 




15 ]/ 20 

(a) Fq: all 3 components (b) F%: all 3 components 

Figure 3. Generic behavior of the rank 3 case . 



(c) F2: all 3 components 



We also find behavior which we conjecture to be parameterized by higher-dimensional 
ultradiscrete hypergeometric functions. For initial conditions and parameters 
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—00 


—00 
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we obtain the behavior exhibited in figure HI 





(a) Fq: all 3 components (b) F\: all 3 components 

Figure 4. Some special behavior of the rank 3 case . 



(c) F 2 : all 3 components 
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6. Conclusions and discussion 

We have presented a noncommutative generalization of g-Prv- Conditions were derived such that 
the matrix valued systems could be ultradiscretized. In section 4, the matrix generalization of 
ultradiscrete Pjy was presented. In section 5, a small snapshot of the rich phenomenology 
was presented. Due to space restrictions, only certain aspects of this phenomenology was 
presented, yet our preliminary findings suggest many avenues for future research, including the 
generalization of the results in [T6j to higher dimensional ultradiscrete hypergeometric functions. 

It is worth noting that a different generalization of q-Pjy, has been studied by Kajiwara et 
al. [12], [13] . It would be interesting to know how both generalizations can be combined. 
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Appendix A. Miscellaneous properties of the group G 

By deducing properties of the group G presented in (13.81) . we may deduce properties of our 
elements {s^ and {zj} since the must be homomorphic images of the generators of G, while 
the {zi} are determined by the {s^ via (13.41) . 

Proposition 1 

So 6 = Si = 9l (A.l) 
Proof Constraint ( 13. 5p implies 

2 — 2 4 4 —2 

9 2 =9i 9 = 9i9 ■ 
Therefore g® = gf, and similarly we have the full proof. □ 
(Note that this implies [g^g^\ = 0.) 

Proposition 2 Group elements {g{\ have order 18. 

Proof As g\ = g® it follows from constraints (13.51) and (13.61) that 



Hence 



8 2 2 -2 

9o = 929o92 ■ 



So 4 = 9l9o92 2 = 9o 



and therefore 

^ 8 = I- (A.2) 

The proofs of 

gl 8 = i , gl 8 = i 

proceed in the same manner. □ 
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Proposition 3 

[gf,g-] = o. (A.3) 

Proof Using (1A.ll) . equation (13.71) shows us that 
9o = 9l 3 9o 3 9l 3 9o 3 , 



further application of (IA.1I) reveals 



00 = 9 e gbo 3 9i 3 , 



hence, using ( 1A.2I) . 



9o 3 9 3 i = 9 3 i9o 3 - 



Therefore we have the commutativity of g$ and gf, and similarly we obtain (1A.3|) . □ 
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